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Symmetry adapted bases in quantum chemistry and bases adapted to quantum informa- 
tion share a common characteristics: both of them are constructed from subspaces of the 
representation space of the group 5*0(3) or its double group (i.e., spinor group) SU{2). 
We exploit this fact for generating spin bases of relevance for quantum systems with 
cyclic symmetry and equally well for quantum information and quantum computation. 
Our approach is based on the use of generalized Pauli matrices arising from a polar de- 
composition of SU{2). This approach leads to a complete solution for the construction 
of mutually unbiased bases in the case where the dimension d of the considered Hilbert 
subspace is a prime number. We also give the starting point for studying the case where 
d is the power of a prime number. A connection of this work with the unitary group U (d) 
and the Pauli group is brielly underlined. 
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Introduction 



The notion of symmetry adapted functions (or vectors) in physical chemistry and sohd 
state physics goes back to the fifties^. The use of bases consisting of such functions 
allows to simplify the calculation of matrix elements of operators and to factorize the 
secular equation. Symmetry adaptation generally requires two type of groups: the sym- 
metry group for the hamiltonian (often a finite group when dealing with molecules) and a 
chain of classification groups for the operators and state vectors (often continuous groups 
like unitary groups^'^'^ and finite groups^'^'^'^'^). The interest of symmetry adapted bases 
(atomic orbitals, molecular orbitals, spin waves, etc.) is well-known in quantum chem- 
istry. In particular, the spherical harmonics (e.g., in atomic spectroscopy) and cubical, 
tetragonal or trigonal harmonics (e.g., in crystal-field theory and ligand field theory^°) are 
quite familiar to the practitioner in theoretical chemistry and chemical physics. 

The symmetry adapted functions generally span bases for finite-dimensional Hilbert 
spaces associated with reducible or irreducible representations of a symmetry group. In 
the case of low dimensions, such spaces are especially useful in the emerging fields of 
quantum information and quantum computation (quantum state tomography and quan- 
tum cryptography), two fields at the crossing of informatics, mathematics and quantum 
physics. In fact, a Hilbert space of finite dimension d can describe a system of qudits 
(qubits correspond to c? = 2, qudits to d arbitrary). Qudits can be realized from many 
physical systems. We undersee that qudits could be also produced from chemical sys- 
tems. 

It is the object of this paper to construct bases which play an important role for quan- 
tum systems with cyclic symmetry and for quantum measurements and quantum informa- 
tion theory. 

The organisation of this paper is as follows. Section 1 is devoted to an alternative 
to the quantization scheme of angular momentum. In Section 2, this scheme 

is worked out for generating bases in a form adapted to physical and chemical cyclic 
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systems as well as to quantum inforaiation. Section 3 deals with somme examples in low 
dimensions. Finally, we develop in Section 4 a systematic construction of generalized 
Pauli matrices which are at the origin of generalized spin bases. In the closing remarks, 
we mention the interest of this work for the special unitary group and the Pauli group. 

Throughout the present work, we use the Dirac notation familiar in quantum chem- 
istry. As usual, stands for the Hermitean conjugate of the operator A. In addition, 
[^4, B]_ and [^4, B]^ denote the commutator and the anticommutator of A and B. Finally, 
i is the pure imaginary. 

1 AN ALTERNATIVE TO THE {fj,} SCHEME 

Let us consider a generalized angular momentum. We note its square and jz its z- 
component. The common eigenvectors of and jz are denoted as \j,m). We know 
that^^ 

m) = j{j + l)|j, m), jz\j, m) = m\j, m) (1) 

in a system of units where the rationalized Planck constant is equal to 1 . For a fixed value 
of the quantum number j (with 2j e N), we note E{2j + 1) the (2j + 1) -dimensional 
Hilbert space spanned by the basis 

= m) : m = j, j - 1, • • • , -j). (2) 

The basis hg is adapted to spherical symmetry (adapted to the group SO{Z) if j is an 
integer or the group SU (2) if j is an half on an odd integer). We take the basis hg in an 
orthonormal form, i.e., the scalar product (j, m\j, m') satisfies 

{j,m\j,m') ^ Sm,m' (3) 

for any value of m and m'. 

In the applications to quantum chemistry, the generalized angular momentum can be 
an angular momentum, a spin angular momentum, a total (spin + orbital) angular mo- 
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mentum, etc. The vectors |j, m) can have several realizations. For instance, in the spec- 
troscopy of 4/^ lanthanide ions, we have state vectors of type | J, M) = \Af^TSLJM) 
in the Russell-Saunders coupling (here j = J and m = M). This constitutes one of many 
possible realizations of the vectors | j, m). 

Besides the basis hs, another interesting basis can be obtained as follows. Let us 
consider the operator 



(4) 



m=-j 



where we use the notation of Dirac for projectors. In Eq. (Hj), we have 

2m 



r G M, a = 0, 1, ■ ■ ■ , 2j, q = exp 



(5) 



2j + 1 

The operator Vra is an extension of the operatorf/^ defined in a previous work^^ (Ur = f ro)- 
From Eq. (H)), we can check that the action of Vra on the state | j, m) is given by 



'ra\j, m) = (1 - 6m,) q^'-'^^'^lj, m + 1) + 5„,,e^2-^-nj, -j). 



Furthermore, the matrix K-a of the operator Vra on the basis bs reads 



(6) 



/ 




Q Q 



\ 



2ja 

/ 



(V) 



where the lines and columns are labeled in the order |j, j), — 1), ■ ■ ■ , \j, — j). It can 
be shown that the operator and Vra commute so that the complete set {j^,fra} of 
commuting operators constitutes an alternative to the set {j'^,jz}- 

We may ask what are the analogues of the vectors \j,m) in the scheme {j^,fra}? 
Indeed, they are the common eigenvectors of the operators and v^a- As a result, these 
eigenvectors are 



\ja; ra) 



y2j + 



1 ^ 

— y 1 



{j+'m){j—'m+l)a/2—jmr+{j+Tn)a 



\j,m), a = 0,l,---,2j (8) 



m=-j 



for a = 0, 1, ■ ■ ■ , 2j. More precisely, we have the eigenvalue equations 

Vra\ja;ra) = g^(''+''^""|ja; ra), f\ja;ra) = j{j + l)\ja;ra). (9) 
For fixed j and a(2j G N, a = 0, 1, ■ ■ ■ , 2j), the basis 

Bra = {\ja-ra): a = 0,1, ■■■,2j} (10) 
is an orthonormal basis since 

{ja;ra\ja';ra) = Sa,a' (H) 

for any value of a and a'. In the particular case where 2j + 1 is a prime integer, the 
overlap between the bases Bra and Brb is such that^^ 

\{ja;ra\jP;rb) \ = 6^^^5a,b + -^^=={1 - Sa,b) (12) 

a property of considerable importance in quantum information. Note that Eq. (fT2l) is 
compatible with Eq. (fTTI) . 

2 A FORMULATION FOR rf-DIMENSIONAL QUAN- 
TUM SYSTEMS 

The parameter r is of interest for group-theoretical analyses but turns out to be of no 
concern here. Therefore, we shall restrict ourselves in the following to the case r = 0. In 
addition, we shall adopt the notation 

k = j — m, \k) = \j,m), \aa) = \ja;Oa), d = 2j + 1 (13) 

that is especially adapted to quantum information (the vectors \0),\1), ■ ■ ■ ,\d—l) are then 
called qudits, the case d = 2 corresponding to ordinary qubits) and to cyclic chemical 
systems (for which \d) = |0), \d+ 1) = . . . ). 
The basis bg becomes 

Ba = {\k):k = 0,l,---,d-l} (14) 
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known as the computational basis in quantum information theory. The action of Vra on 
the basis Ba of S{2j + 1) is described by 

Voa\k) = q''^\k-l) (15) 

where k — 1 should be understood modulo d (i.e., \ — 1) = \d — 1)). The vectors \aa) of 
the orthonormal basis 

Bra = {\aa) : a = 0,1, ■ ■ ■ ,d — 1} (16) 

can be written as 

\aa) = ^ Vg('^'^')('+')"/'-('=+')"|A;) (17) 

where a can take the values a = 0,1, - ■ ■ ,d — 1. These vectors satisfy the eigenvalue 
equation 

wojaa) = g('^-^)"/2-"|aa) (18) 

that corresponds non a nondegenerate spectrum for the operator foa- 

All relations given in Section 1 up to this point are valid for d arbitrary. In the special 
case where dis a prime integer, Eq. ([T2)) yields 

\{aa\b(3)\ = 5^,f,5a,t + ^(1 - Sa,b) (19) 

a relation valid for any value of a, b, a and /3 in the set {0, 1, ■ ■ ■ , c? — 1}. In quantum in- 
formation, two bases Boa and Bqi, satisfying Eq. (fT9l ) are said to be mutually unbiased^^. 
Such bases play an important role in quantum cryptography and quantum state tomogra- 
phy. It is well-known that a complete set ofd+l mutually unbiased bases can be found 
when is a prime integer or the power of a prime integer. 
We continue with some typical examples. 
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3 SOME TYPICAL EXAMPLES 



The case d = 2 

In this case, relevant for a spin j = 1/2 or for a qubit, we have q = —1 and a,aE {0, 1}. 
The matrices of the operators Voa are 

We note in passing a connection (to be generalized below) with the Pauli matrices since 
Vbo = and Vqi = —iuy. From Eqs. (fT4l) . (fT6l ) and (fTTI) . the bases i?2, -Boo and i?oi are 

B2 : |0), |1) (21) 
^00 : |00) = i=(|0) + |l)), |01) = i=(-|0) + |l)) (22) 

5oi : |10) = i=(z|0) + |l)), |ll) = -L(-z|0) + |l)) (23) 

which satisfy Eq. (fT9l ). Note that by using the spinorbital 

o^=\\.\) = \^). /5=l^,-^) = |l) (24) 

(a for spin up and (3 for spin down) familiar to the quantum chemist, Eqs. (l2T1)- (|23l) can 
be rewritten as 

B2 : a, (3 (25) 

Boo : |00) = -^(a + /3), |01) = --^ (a - /3) (26) 

Boi : |10)=z-^(a-z/5), |11) = -z-^ (a + z/3) . (27) 

In terms of eigenvectors of the matrices Voa, we must replace the vectors \aa) by column 
vectors. This leads to 



: 100) 101) --^(_\) (29) 
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The case d = 3 

This case corresponds to a spin j = 1 or to a qutrit. Here, we have q = exp{2ni/3) and 
a, a G {0, 1, 2}. The matrices of the operators voa are 



Vr 



00 



1 

1 

1 



q 

g2 

1 



g2 

g 

1 



(31) 



The bases B^, Bqq and Bqi B02 are 



B3 
Boo 



01 



5, 



02 



0), |1), |2) 

00) = ^ (|0) + |1) + |2)) , |01) = (g^|0) + g|l) + |2)) 



02) 



(g|0) + g^|l) + |2)) 



10) = ^ (g|0) + g|l) + |2)) , 111) = ^ (|0) + q'\l) + |2)) 
12) = ^(g^|0) + |l) + |2)) 

20) = ^ (g^|0) + q'\l) + |2)) , |21) = i= (g|0) + |1) + |2)) 
22) = -^(|0) + g|l) + |2)). 



They satisfy Eq. (fT9l ). In terms of colum vectors, we have 



5, 



00 



5, 



01 



5, 



02 




(32) 
(33) 

(34) 

(35) 

(36) 

(37) 

(38) 



(39) 
(40) 
(41) 
(42) 
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The case (i = 4 

This case corresponds to a spin j = 3/2. Here, we have q = i and a, a G {0, 1, 2, 3}. 
Equation (flTI) can be applied to this case too. However, the resulting bases B4,, Bqq, Bqi, 
Bq2 and i?03 do not constitute a complete system of mutually unbiased bases (rf = 4 is 
not a prime number). Nevertheless, it is possible to find d + I = 5 mutually unbiased 
bases because = 2^ is the power of a prime number. This can be achieved by replacing 
the space £^(4) spanned by {|3/2,m) : m = 3/2, 1/2, —1/2, —3/2} by the tensor product 
space £(2) £(2) spanned by the basis 

{a® (43) 

The space £{2) £(2) is associated with the coupling of two spin angular momenta 
ji = 1/2 and = 1/2 or two qubits (in the vector v, u and v correspond to ji and 
j2, respectively). An alternative basis for £{2) ® £{2) is 

{a® a,]^{a® (3 + (3 ® a),(3 ® (3,]^{a® (3 - (3 ® a)}. (44) 

The vectors in (|43l) are well-known in the treatment of spin systems. The first three vectors 
are symmetric under the interchange 1^2 and describe a total angular momentum J = 1 
while the last one is antisymmetric and corresponds to J = 0. It should be observed 
that the basis (|44|) illustrates a connection between the special unitary group SU (2) and 
the permutation group S2 (a particular case of a reciprocity theorem between irreducible 
representation classes of SUn and Sm)- 

In addition to the bases (l43l) and (|44)) . it is possible to find other bases of £{2) ® £{2) 
which are mutually unbiased. The d = A mutually unbiased bases besides the canonical 
or computational basis (|43] ) can be constructed from the eigenvectors 

\abal3) = \aa) (45) 

of the operators Wab = ® ^"06 (the vectors \aa) and \bf3) refer to the two spaces £{2)). 
As a result, we have the d+l = 5 following mutually unbiased bases where A = {l — i)/2 
and /i = {l+i)/2. 
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The canonical basis: 

a® a, a<S)/3, (3 <S) a, (5 ® (5 

or in column vectors 

/1\ /o\ /0\ /0\ 

10 

' ' 1 ' ■ 

VI VI VI VI 

The wqq basis: 

|0000) = ]-{a®a + a®l3 + l3®a + l3®l3) 

|0001) = + 

|0010) = \-{a®a^a®i5-i5®a-i5®i5) 

|0011) = \-{(x®a- a® ii - i5®a\ 

Zi 

or in column vectors 









i ^ \ 








/ 1 \ 


1 


1 


1 


-1 


1 


1 


1 




2 


1 


' 2 


1 


' 2 


-1 


^ 2 


-1 




V) 




\-^) 




\-'^) 




^ 1 / 



The wii basis: 

1 1100) = ^{a(S)a + ia®) P + if3 ®a - f3 ® f3) 

|1101) = ^{a(S)a-ia(S) P + ip<S>a + p® 
z 

|1110) = ha®a + ia® P -ip®a + P® P) 
z 

jllll) = ha<S)a-ia®)P-iP<S)a- P<S)P) 
z 

or in column vectors 





( ^ ^ 




( ^ \ 




/ 1 ^ 




/ 1 \ 


1 


i 


1 


—i 


1 


i 


1 


—i 


2 


i 


' 2 


i 


' 2 


—i 


' 2 


—i 




V^l 




\ 1 ) 




\ 1 ^ 




V-iy 
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A|0100) + /x|0111) = 




$ q; + q; $ 






5/3) 


(58) 


/x|0100) + A|0111) = 


I, 


^ a — 


^l3 + il3(i 


^a + ip^ 


^P) 


(59) 


A|0101) +/x|0110) = 


I, 


^ a — 




^ a — ip <i 


^P) 


(60) 


/x|0101) + A|0110) = 


1, 

2(«^ 


^ a + 


^f3 + ip^ 


^ a — ip <^ 


^P) 


(61) 



or in column vectors 





( 1 ^ 




/ 1 \ 




( ' ) 




( 1 \ 


1 


1 


1 


-1 


1 


-1 


1 


1 


2 


—i 


' 2 


i 


' 2 


—i 


' 2 


i 




\ ^ J 




\ ^ J 




\-0 







The wio basis: 



or in column vectors 

/ 1 \ 



1 



/ 1 \ 

i 

-1 

V ^ / 



/ 1 \ 

i 
1 



/ 1 \ 

—i 
-1 



(62) 



A| 1000) 1011) = 


1, 


^ a — ia<i 


+ 


^a + iP^ 


5/3) 


(63) 


//|1000) + A|1011) = 


1, 


^ a + ia i 


^P-P^ 


^a + ip^ 


5/3) 


(64) 


A|1001) + /i|1010) = 


1, 


§ a + ia 


^P + P<i 


^ a — ip ^ 


5/3) 


(65) 


/x|1001) + A|1010) = 


1, 


$ a — ia 


^P-P^ 


^ a — ip Q 


5/3) 


(66) 



(67) 



It is to be noted that the vectors of the woq and wn bases are not intricated (i.e., each 
vector is the direct product of two vectors) while the vectors of the wqi and wio bases are 
intricated (i.e., each vector is not the direct product of two vectors). 

4 GENERALIZED PAULI MATRICES 



From the operators voa, it is possible to define two basic operators x and z which can be 
used for generating generalized Pauli matrices. Let us put 



X = VoQ, z = vIqVoi. 
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(68) 



The action of x and z on the space £{2j + 1) is given by 

x\j, m) = (1 - 6^j) I J, m + 1) + 6m,j\3, -j) ^ x\k) = \k - 1) (69) 

and 

z\j, m) = q^-'^lj, m) ^ z\k) = q''\k) (70) 

where q = exp{27ri/d) with d = 2j + 1. The ci-dimensional matrices X and Z of x and 
z are 



x= ,z= oog^--- o . (vi) 



/o 


1 


■ 


• ^\ 




/I 





■ 


■ ° \ 








1 ■ 


■ 







Q 


■ 













, z = 








q' ■ 











• 


• 1 












\1 





• 






Vi 





■ 





The operators x and z are unitary and satisfy the g-commutation relation 

xz — qzx = 0. (72) 

Pairs of operators satisfying a relation of type (l72l) are nowadays refered to as Weyl pairs. 
Thus, the operators x and z constitute a Weyl pair. Weyl pairs were introduced at the 
begining of quantum mechanics^^. They were extensively used for factorizing the secular 
equation in connection with a study of alternating hydrocarbons^^ and for constructing 
analogues of the usual Pauli matrices^^'^'''^^'^^. 
Let us now define the operators 



Uab = x°'z'^, a,b = 0,1, ■ ■ ■ , d — 1. 



The d'^ operators Uab are unitary and satisfy the following trace relation 



(73) 



Tr£-(2j+l) (ulf^Ua'b'^ = (2j + 1) Sa,a' Sb,b' 



(74) 



where the trace is taken on the rf-dimensional space £^(2j + 1). Additionally, the commu- 
tator [uab, Ua'b']- and the anti-commutator [uab, Ua'b'] + of Uab and Ua'b' are given by 



[Uab, Ua'b'\:f 



l^q-ba' ^ q-ab'^ ^ a" = a + a', b" = b + b'. 



(75) 
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Consequently, [uab, Ua'b']- = if and only if ab' — ba' = (mod d) and [uab, Ua'b']+ = 
if and only if ab' — ba' = {l/2)d (mod d). Therefore, all anti-commutators [uab,Ua'b'] + 
are different from if d is an odd integer. 

Two consequences follow from Eqs. (1741 ) and (|75l ). First, the trace relation (1741 ) shows 
that the rf^ operators Uab are pairwise orthogonal operators so that they can serve as a basis 
for developing any operator acting on the Hilbert space £{d). Second, the commutation 
relation (fTSl) shows that the set {uab : a,b = 0,1, ■ ■ ■ , d — 1} generates a rf^-dimensional 
Lie algebra. This algebra turns out to be the Lie algebra of the unitary group U{d). The 
subset {uab '■ a,b = 0,1, - ■ ■ ,d — 1} \ {mqo} thus spans the Lie algebra of the special 
unitary group SU{d). 

All this is reminiscent of the group SU{2), the generators of which are the well-known 
Pauli matrices. Therefore, the operators Uab shall be refered as generalized Pauli operators 
and their matrices as generalized Pauli matrices. As an illustration, let us deal with the 
cases d = 2 and d = 3. 

Exemple 1 

In the case j = 1/2 ^ d = 2 (^ q = —1), the matrices of the 4 operators Uab with 

a, 6 = 0, 1 are 

/ = X°Z°=(J X = X'Z'=(^l J) (76) 

Z = X'Z'=(l -l)' Y = X'Z'=(l "J). (77) 

In terms of the usual (Hermitean and unitary) Pauli matrices a^, ay and cr^, we have 
X = cr^, F = —iay, Z = Oz- The approach developed in the present paper lead to Pauli 
matrices in dimension 2 that differ from the usual Pauli matrices. This is the price one has 
to pay in order to get a systematic generalization of Pauli matrices in arbitrary dimension. 
It should be observed that the commutation and anti-commutation relations given by (1751) 
with d = 2 correspond to the well-known commutation and anti-commutation relations 
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for the usual Pauli matrices (transcribed in the normalization X^Z^ — ax, X^Z^ — —iuy, 

Exemple 2 

In the case j — 1 d — ?> q — exp(27ri/3)), the matrices of the 9 operators Uab with 
a,b = 0, 1, 2, viz., 

= I X^Z^ = X X^Z^ = X'^ X^Z^ = Z X^Z'^ = (78) 



X^Z^ = XZ X'^Z^ X^Z^ = X^Z X^Z^ = XZ'^ (79) 



are 





X= 1, X2= l (80) 




1 0\ (\ ^ 0' 

Z=|0 g , Z^^ O g2 o|, XZ^ \ g^l (81) 
?V \o 5 




^2^2 _ I ^ Q Q I jsc^^ = I 1 1, = I q \ . (82) 






These generalized Pauli matrices differ from the Gell-Mann matrices^^ used in elementary 
particle physics. They constitute a natural extension of the Pauli matrices in dimension 

5 CONCLUDING REMARKS 

The various bases described in the present paper are of central importance in quantum 
information and quantum computation. They also play an important role for quantum 
(chemical and physical) systems with cyclic symmetry. By way of illustration, we would 
like to mention two examples. 
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Let us consider a ring shape molecule with N atoms (or agregates) at the vertices of 
a regular polygon with sides (N = 6 for the benzen molecule CqRg). The atoms are 
labelled by the integer n with n = 0,l,---,A^ — 1. Hence, the cyclic character of the ring 
shape molecule makes it possible to identify the atom with the number n to the one with 
the number n + kN where k E Z (the location of an atom is defined modulo N). Let \Lpn) 
be the atomic state vector, or atomic orbital in quantum chemistry parlance, describing 
a 7r-electron located in the neighboring of site n. From symmetry considerations, the 
molecular state vector, or molecular orbital, for the molecule reads^^ 



with s = 0, 1,---,A^— 1. Asa result, the molecular orbital l^^) assumes the same form, 
up to a global phase factor, as the state \aa) given by Eq. (flTI) with a = and a = s. 

A similar result can be obtained for a one-dimensional chain of N 1/2-spins (num- 
bered with n = 0, 1, ■ • ■ , — 1) used as a modeling tool of a ferromagnetic substance. 
Here again, we have a cyclical symmetry since the spins numbered n = N and = are 
considered to be identical. The spin waves can then be described by state vectors^^ very 
similar to the ones given by Eq. (flTI) with again a = 0. 

We close this work with two remarks of a group-theoretical nature, one concerning a 
continuous group, the other a finite group, connected with the operators Uab- 

First, as mentioned in Section 4, the set {uab : a,b = 0,1, ■ ■ ■ , d—l}\{uoo} constitutes 
a basis for the Lie algebra SU{d). Such a basis differs from the well-known Cartan basis 
or from the Gel'fand-Tsetlin basis. In the special case d = p, with p prime integer, the 
basis {uab : a,b = 0,1, ■ ■ ■ , d — 1} \ {uqq} can be partioned into p + 1 disjoint subsets, 
each subset containing p — 1 commuting operators^^'^^. In other words, it is possible to 
decompose the Lie algebra of SU (p) into p + 1 Cartan subalgebras of dimension p — 1. 
It can be proved that each subalgebra is associated with a basis of £(p) and that the set 
of the p + 1 corresponding bases is a complete set of mutually unbiased bases. A similar 
decomposition holds of SU (d) in the case where d = p^, with p prime integer and e 




N-l 



(83) 



n=0 
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positive integer^^. However, in this case we need to replae S{d) by £{p)®^. 

A second group-theoretical remark concern a finite group known as the Pauli group 
or the finite Heisenberg-Weyl group^^'^®'^^'^^'^^. The set {uab : a, 6 = 0, 1, • • • , ci - 1} is 
not closed under multiplication. However, it is possible to extend the latter set in order to 
have a group. For this purpose, let us define the operators Wabc via^^ 

Wabc = (fuhc, a,6, c = 0, 1. (84) 

Then, the set {wahc : a, & = 0, 1, • • • , — 1}, endowed with the multiplication of opera- 
tors, is a group of order cP. This group (the Pauli group) is of paramount importance in 
quantum information and quantum computation^^'^^. 
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